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Executive  Summary 


A  computationally  fast,  mesoscale  meteorological  analysis  and  very  short  range 
prediction  (nowcasting)  capability  has  been  developed  for  the  U.S.  Army,  specifically  for 
the  Integrated  Meteorological  System  (IMETS),  which  is  the  Army’s  tactical  weather 
system.  This  tactical  analysis  and  nowcast  capability  is  currently  referred  to  as  the 
Weather  Running  Estimate-Nowcast  (WRE-N)  by  the  U.S.  Anny  Research  Laboratory, 
White  Sands  Missile  Range,  NM,  which  is  serving  as  the  Army  lead  in  its  development. 
The  WRE-N  discussed  in  this  report  is  an  initial  implementation  of  such  a  capability, 
configured  for  a  domain  of  approximately  100  km  by  100  km,  with  a  horizontal  grid 
resolution  of  1.5  km.  Additionally,  this  WRE-N  has  been  tailored  to  accept 
meteorological  data  sources  routinely  available  to  the  IMETS,  including  Air  Force 
Weather  Agency  (AFWA)  Mesoscale  Model  Version  5  (MM5)  numerical  forecasts, 
surface  sensor  observations,  and  upper-air  radiosonde  observations.  Finally,  this 
implementation  allows  three-dimensional  gridded  distributions  of  temperature,  dew-point 
temperature,  and  horizontal  wind  vector  components  to  be  analyzed  and  nowcasted. 

The  method  is  evaluated  over  two  different  geographic  areas:  one  centered  near  Dallas, 
TX,  and  the  other  near  Denver,  CO.  This  report  shows  that  the  method  appeared  to  add 
value  (when  compared  to  AFWA  15-km  MM5  background  forecasts)  for  up  to  3  hr  over 
the  Dallas  domain  (fairly  uniform  terrain)  and  up  to  2  hr  over  the  Denver  domain 
(complex  terrain).  Although  further  studies  are  needed  over  different  areas  and  seasons, 
we  tentatively  conclude  that  the  present  WRE-N  method,  which  is  based  on 
extrapolation,  can  be  used  for  at  least  two  hours  for  many  meteorological  scenarios, 
especially  when  the  background  forecast  model  (the  operational  AFWA  MM5)  is  of  good 
quality  on  the  mesoscale. 


Intentionally  Left  Blank. 


1.  Introduction 


For  U.S.  Army  operations,  being  able  to  accurately  update  high  spatial  and  temporal 
resolution  meteorological  information,  generated  by  fusing  or  assimilating  current  local 
meteorological  observational  data,  is  more  desirable  than  relying  on  results  obtained  by 
previously  executed  model  forecasts  alone.  Although  current  state-of-the-art  numerical 
weather  prediction  (NWP)  models  offer  an  unprecedented  capability  to  simulate  small- 
scale  mean  atmospheric  structures,  individual  detenninistic  model  solutions  tend  to  suffer 
from  day-to-day  spatial  and  temporal  errors  in  amplitude,  phase,  and  variability.  These 
errors  are  due  to  many  reasons,  whose  discussion  is  beyond  the  scope  of  this  report. 

To  correct  for  these  kinds  of  errors  and  model  biases,  it  is  necessary  to  assimilate  real 
meteorological  observational  data  through,  for  example,  NWP  model  forecast  cycling. 

The  advanced  assimilation  methods  that  have  been  developed  over  the  past  few  decades 
now  offer  dynamically,  physically  consistent  and  statistically  optimal  methods  of 
“correcting”  NWP  model  tendencies  on  the  fly.  If  the  “current”  weather  situation  is 
accurately  analyzed,  it  may  be  possible  to  reasonably  predict  very  short-term  weather 
situations  (nowcasting).  Our  ability  to  accomplish  this  will  increase  as  techniques 
continue  to  improve  and  better  methods  of  assimilating  observed  cloud-scale  fields  and 
dynamically  coupling  them  to  NWP  model  fields  are  produced. 

Rapid  adjustments  of  intelligence  and  warfighter  decision  making  tools,  to  account  for 
changes  in  the  local  battlefield’s  three-dimensional  (3D)  environmental  state  (which  often 
deviate  from  the  existing  operational  NWP  model  forecast  guidance),  are  known  as 
Weather  Running  Estimates  (WREs).  For  the  Anny,  near-surface  meteorological  changes 
on  timescales  ranging  from  minutes  to  a  few  hours  are  critical  to  ground  operations.  Such 
changes  are  classified  as  meso-y  (1-  to  10-km  grid  spacing)  scale  phenomena  (Orlanski, 
1975).  This  range  is  also  commonly  referred  to  as  the  “cloud”  or  “storm”  scale. 

Nowcasting  is  a  method  of  diagnosing  the  current  weather  situation  by  combining 
existing  forecast  guidance  with  observational  data,  using  various  methods  of 
analysis/fusion  and  assimilation,  and  then  extending  this  meteorological  “picture” 
forward  in  time  to  produce  a  very  short-range  forecast  or  prediction  (Browning,  1982). 
The  predictive  component  can  be  produced  in  many  ways,  such  as  through  cycling 
(coupling  analyses  to  prognostic  NWP  models),  extrapolation,  statistical  methods,  or 
expert  system  approaches.  The  Weather  Running  Estimate-Nowcast  (WRE-N) 
incorporates  a  3D  gridded  analysis  (the  WRE)  and  a  3D  short-range  gridded  prediction 
counterpart  (the  Nowcast). 

Using  conventional  observational  data  along  with  operational  Mesoscale  Model  Version 
5  (MM5)  NWP  model  data  produced  at  the  U.S.  Air  Force  Weather  Agency  (AFWA),  a 
relatively  simple  and  fast  (in  terms  of  execution)  WRE-N  system  has  been  developed, 
building  upon  earlier  analysis  techniques  constructed  at  the  U.S.  Army  Research 
Laboratory  (ARL),  White  Sands  Missile  Range,  NM  (Henmi,  2003).  In  the  current  WRE- 
N  system,  separate  surface  and  upper-air  meteorological  analyses  are  first  produced 
univariately.  To  produce  a  surface  analysis,  local  surface  observations  and  concurrent 
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NWP  model  forecast  data  are  composed  using  a  successive  correction  method  (Sashegyi 
and  Madala,  1994).  Similarly,  the  upper-air  analyses  (produced  for  many  predetermined 
levels  above  ground)  are  a  composition  of  upper-air  radiosonde  soundings  and  NWP 
model  forecast  data,  using  the  same  scheme.  After  both  the  surface  and  upper-air 
analyses  have  executed,  similarity  theory  formulae  named  MARIAH  (Rachele  et  al., 
1995)  are  used  to  couple  the  surface  and  upper-air  analyses  into  a  single  set  of  3D  fields. 
This  3D  “snapshot”  is  the  WRE  component  of  the  WRE-N  system.  The  analysis  method 
used  to  produce  the  WRE  has  been  previously  statistically  evaluated  (Hemni  and  Dumais, 
2004)  over  Utah  using  the  University  of  Utah’s  MesoWest  data  (Horel  et  al.,  2002).  That 
study  showed  that  this  method  produced  much  better  statistical  results  for  temperature, 
dew-point  temperature,  and  wind  vectors  as  compared  to  previously  generated 
background  forecast  fields  from  the  Pennsylvania  State  University/National  Center  for 
Atmospheric  Research  MM5. 

The  U.S.  Army’s  Integrated  Meteorological  System  (IMETS)  regularly  receives  (from 
AFWA)  surface  observation  data  (such  as  METAR),  World  Meteorological  Organization 
(WMO)  upper-air  radiosonde  sounding  data,  and  WMO  gridded  binary  (GRIB)  data 
produced  from  operational  forecasts  of  the  MM5.  The  MM5  forecast  data  are  typically 
produced  twice  daily  at  a  horizontal  grid  spacing  of  15  km.  Using  the  AFWA  operational 
MM5  fields  for  both  current  and  “future”  background  field  data,  the  WRE-N  method  has 
been  developed  and  adapted  for  the  IMETS. 

To  produce  the  Nowcast  component  of  the  WRE-N,  the  successive  corrections  analysis 
method  has  been  extended  to  make  a  short-range  “prediction”  for  the  period  out  to  3  h 
forward.  The  WRE-N  has  been  nominally  configured  for  a  domain  of  approximately 
100  km  by  100  km  and  for  a  horizontal  grid  spacing  of  about  1.5  km.  For  clarification, 
the  successive  corrections  approach  is  used  to  produce  the  WRE  “current”  analysis,  while 
a  linear-interpolation  method  is  used  to  extrapolate  the  WRE  fields  forward  to  produce 
the  short-range  nowcast.  In  extrapolating  the  WRE  analysis  forward,  the  AFWA  MM5 
fields  (valid  at  the  desired  forward  times)  are  used  as  background  data.  The  method  has 
recently  been  applied  to  two  different  geographical  domains  and  evaluated  statistically, 
with  the  preliminary  results  presented  later  in  this  report. 


2.  Data  Used 


The  WRE-N  is  designed  to  incorporate  the  data  most  readily  available  in  the  current 
IMETS  system,  which  includes  digital  terrain  elevation  data  and  both  observational  and 
NWP  model  forecast  data  supplied  from  AFWA  (surface,  upper-air  radiosonde,  and 
GRIB-formatted  MM5  output). 
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2.1  Elevation  Data 

The  elevation  data  file,  which  encompasses  the  WRE-N  area  of  interest  (AOI),  is 
produced  by  reading  the  Digital  Terrain  Elevation  Data  (DTED)  Level  0  from  the 
Defense  Mapping  Agency  (DMA).  An  existing  IMETS  computer  program  produces  a 
terrain  elevation  data  file  for  a  161  -by- 161  gridded  domain  centered  at  the  desired 
WRE-N  AOI,  with  the  data  interpolated  from  the  native  30  arc  second  resolution  to  the 
desired  WRE-N  1.5-km  grid  spacing.  This  is  done  on  a  pseudo-Universal  Transversal 
Mercator  (UTM)  grid,  which  is  the  same  grid  used  by  the  WRE-N.  Since  the  generated 
elevation  data  covers  a  domain  somewhat  larger  than  that  required  for  the  WRE-N 
(100  km  by  100  km),  the  actual  terrain  data  used  by  WRE-N  is  clipped  from  the  original 
file  of  16 l-by-16 1  grid  points. 

2.2  MM5  Data 

From  the  GRIB-formatted  MM5  forecast  data  files  sent  from  the  AFWA,  vertical 
distributions,  or  “model  soundings,”  of  forecast  fields  at  3-h  output  intervals  are  extracted 
for  the  AOI  region  and  placed  into  ASCII  input  files.  These  files  include  MM5  data  that 
are  interpolated  to  both  the  desired  WRE-N  “current”  hour  and  3  hr  forward  from  that 
time.  The  base  time  of  the  MM5  must  be  between  9  and  24  hr  earlier  than  the  WRE 
0  hour.  A  portion  of  the  data  contained  within  the  WRE-N  input  MM5  ASCII  files  is 
shown  in  figure  1 . 
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Figure  1.  A  portion  of  an  ASCII  MM5  data  file  used  as  WRE-N  input,  extracted  from  the  MM5  GRIB 
datasets  transmitted  from  AFWA. 
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The  top  lines  of  each  data  file  contain  the  valid  date  and  forecast  time  of  the  MM5 
forecast  data,  the  MM5  model  grid  resolution,  the  number  of  MM5  grid  point 
“soundings”  in  the  file,  and  the  number  of  MM5  vertical  layers.  Subsequent  lines  include 
the  latitude,  longitude,  and  model  terrain  elevation  for  each  MM5  grid  point  sounding, 
followed  by  the  forecast  data  for  the  meteorological  fields  at  each  model  sigma  level.  The 
data  columns  contain  pressures  (mb),  heights  above  sea  level  (m),  temperatures  (°K), 
dew-point  depressions  (°C),  and  wind  vector  components,  it  and  v  (m/s). 

2.3  Surface  Data 

Surface  meteorological  observations  are  input  into  the  WRE-N  as  ASCII  formatted  files 
and  contain  pressure  (mb),  temperature  (°C),  dew-point  temperature  (°C),  wind  speed 
(m/s),  and  wind  direction  (deg)  data.  Figure  2  shows  a  portion  of  the  WRE-N  input 
surface  meteorological  data  file. 

Pressure  (mb) 

7 

32.680000  -96.860001  22Apr  11:53  1011.90 

32.810001  -97.360001  22Apr  11:53  1012.10 

33.180000  -96.580002  22Apr  11:53  1013.50 

33.200001  -97.180000  22Apr  11:53  1013.10 

Temperature  (C) 

12 

32.560001  -97.300003  22Apr  12:05  17.00 

32.660000  -97.099998  22Apr  12:09  17.00 

33.200001  -97.180000  22Apr  11:53  15.00 

33.650002  -97.199997  22Apr  12:05  12.00 

Dewpoint  Temperature  (C) 

12 

32.560001  -97.300003  22Apr  12:05  16.00 

32.660000  -97.099998  22Apr  12:09  16.00 

33.200001  -97.180000  22Apr  11:53  11.00 

33.650002  -97.199997  22Apr  12:05  9.00 

Wind  Speed  (m/s) 

12 

32.560001  -97.300003  22Apr  12:05  1.54 

32.660000  -97.099998  22Apr  12:09  1.54 

33.200001  -97.180000  22Apr  11:53  3.60 

33.650002  -97.199997  22Apr  12:05  2.06 

Wind  Direction  (deg) 

12 

32.560001  -97.300003  22Apr  12:05  250.00 

32.660000  -97.099998  22Apr  12:09  240.00 

33.200001  -97.180000  22Apr  11:53  .00 

33.650002  -97.199997  22Apr  12:05  340.00 

Figure  2.  A  portion  of  an  ASCII  surface  meteorological  data  file  used  as  WRE-N  input,  containing  pressure, 

temperature,  dew-point  temperature,  wind  speed,  and  wind  direction  for  different  observation  locations. 
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2.4  Upper-Air  Radiosonde  Data 

The  upper-air  radiosonde  observation  data  are  ingested  into  the  WRE-N  as  ASCII 
formatted  files,  similar  to  the  MM5  and  surface  observations.  An  abbreviated  portion  of  a 
WRE-N  input  upper-air  radiosonde  data  file  is  shown  in  figure  3.  Information  about 
latitude;  longitude;  height  above  ground;  date  and  hour  of  observation;  and  number  of 
levels  are  given  in  the  first  three  lines.  Pressure  (mb),  height  above  sea  level  (m), 
temperature  (°C),  dew-point  temperature  (°C),  wind  speed  (kn),  and  wind  direction  (deg) 
data  are  also  given.  All  missing  data  are  expressed  as  “-999. 0.” 


32.830002  -97.300003  196.000000 
20 Apr  00:00 
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Figure  3.  A  portion  of  an  ASCII  upper-air  radiosonde  sounding  data  file  used  as  WRE-N  input,  containing 

pressure,  height  above  sea  level,  temperature,  dew-point  temperature,  wind  speed,  and  wind  direction. 

Figure  4  shows  an  example  of  one  of  the  WRE-N  domains  evaluated  in  this  report, 
depicting  the  terrain  and  meteorological  data  sources  used  as  inputs.  The  elevation 
contours  shown  are  for  a  75-by-75  grid  point  domain  with  a  1.5-km  horizontal  resolution, 
centered  at  32.83°  N.  and  97.30°  W.  The  locations  of  the  MM5  grid  points  (15-km 
resolution)  are  marked  with  an  asterisk  (*),  and  surface  and  upper-air  data  locations  are 
represented  with  an  “S”  and  a  “U,”  respectively.  The  city  of  Dallas,  TX,  is  located  near 
the  center  area  of  this  domain,  referred  to  as  the  Dallas  domain  in  this  report. 


5 


Figure  4.  Elevation  contour  distribution  for  the  domain  of  75-by-75  grid  points  with  1.5-km  grid 
spacing,  centered  at  32.83°  N.  and  97.30°  W.  (Dallas  domain). 

NOTE:  The  symbol  represents  MM5  grid  points,  while  “S”  and  “U”  represent  locations  of 
surface  and  upper-air  radiosonde  meteorological  observations,  respectively. 


3.  Analysis  Methods  to  Produce  WRE 


3.1  3D  Analysis  of  MM5  and  Upper-Air  Radiosonde  Data 

The  15 -km  resolution  MM5  forecast  data,  as  shown  in  figure  1,  are  first  interpolated 
spatially  1)  onto  predetermined  height  levels  above  mean  sea  level  and  2)  to  the  WRE-N 
grid  point  locations  at  the  finer  horizontal  grid  resolution  (1.5  km). 

There  are  differences  between  the  terrain  heights  of  grid  points  interpolated  from  MM5 
forecast  data  and  those  obtained  from  DTED  Level  0.  Three-dimensional  temperature  and 
pressure  field  errors,  due  to  the  terrain  height  difference,  are  corrected  as  follows: 
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(1) 


Surface  temperature  and  pressure  are  recalculated  as 

JrTl 

T  =7\  +  (— )*az 

1  1  s0  dz 

dP 

P  =  Pft  +  (— )*AZ 
r  dz 

where 

Ts:  corrected  surface  temperature, 

Ps:  corrected  surface  pressure, 

Ts(}:  surface  temperature  interpolated  from  MM5  data, 

Ps0:  surface  pressure  interpolated  from  MM5  data, 
dT/dr.  temperature  gradient  (6.5  °C  per  1,000  m), 
dP/dz :  pressure  gradient  (9.0  mb  per  100  m), 

A Z:  terrain  height  difference,  Zso  -  Zs, 

Zs0:  terrain  height  interpolated  from  MM5  data,  and 
Zs:  terrain  height  from  DTED  Level  0. 

Vertical  temperature  profile  is  recalculated  as 


T{Z)  =  T+{ 


PZ)-T, 

z -zs 


)  *  (z  -  Z)  ■ 


Here  Zc  is  the  height  determined  empirically  to  correct  vertical  temperature  profile. 
Vertical  pressure  profile  is  recalculated  using  the  hydrostatic  equation: 


^(Z2)  =  ^(Z,)*exP 


■(4xz2 


■z,x- 


■  +  ■ 


nz)  pz2) 


) 


(2) 


(3) 


(4) 


After  these  steps  are  accomplished,  the  data  are  then  vertically  interpolated  once  again  to 
the  vertical  coordinate  levels  of  the  WRE-N.  This  sigma-z  vertical  coordinate  of  the 
WRE-N  is  defined  as 


where 

z  =  the  Cartesian  vertical  coordinate, 

2  =  the  ground  elevation, 

H  =  the  material  surface  top  of  the  model  in  the  z*  coordinate,  and 
H=  the  corresponding  height  in  the  z  coordinate  defined  by  H  =  H  +  zgmax 


(5) 
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where 


=  the  maximum  value  of  the  terrain  elevation  in  the  domain  of  the  WRE-N. 

gmax 


Currently,  the  vertical  extent  of  the  model  ( H  )  is  7,000  m  above  the  highest  grid  point 
elevation  of  the  domain. 


The  following  procedures  are  taken  to  interpolate  the  15-km  resolution  MM5  forecast 
data  to  the  3D  spatial  locations  of  the  WRE-N  AOI: 


1 .  MM5  parameters  (temperature,  dew-point  temperature,  and  horizontal  u  and  v  wind 
components)  are  vertically  interpolated  to  30  predetermined  vertical  levels 

(z/,  i=  1,  30). 

2.  Horizontal  interpolation  of  MM5  parameters  (on  z,  levels)  to  WRE-N  grid  point 
locations  are  perfonned  for  each  parameter,  using  a  1/r"  weighting  factor,  where  “r” 
is  the  distance  between  the  MM5  grid  point  and  the  WRE-N  grid  point. 

3 .  An  arbitrary  parameter  cp  on  a  z,  level  is  linearly  interpolated  to  a  z  *  level  of  the 
WRE-N  grid: 

<p^=<p^P^)L-c)  <«) 

*+i  ~>  k' 


where  ^  is  the  Cartesian  height  above  sea  level  of  the  z  level,  calculated  from 
equation  1  as 


:Ze  +  Z 


t H 


+z 


gmax 


H 


(7) 


and  is  the  height  of  the  kth  z,  level. 


If  there  are  no  upper-air  radiosonde  sounding  data  available  for  a  desired  time,  the  3D 

* 

data  of  cp(z  )  (as  detennined  by  equation  6)  are  used  as  the  3D  background  dataset  for 

the  creation  of  the  WRE-N.  For  the  nowcast  periods  of  the  WRE-N  (to+30,  to+60,  to+90, 
to+120,  /o+l  50,  to+180  min),  this  is  how  background  MM5  fields  are  always  generated, 
since  no  radiosonde  data  are  available  at  these  times.  If  there  are  local  upper-air 
radiosonde  sounding  data  available  at  the  WRE/analysis  time  (to),  the  3D  MM5  data 
calculated  in  equation  2  are  also  composed  with  observed  upper-air  radiosonde  sounding 
data  as  follows: 


1 .  At  each  radiosonde  location,  sounding  data  (such  as  horizontal  wind  vector 
components,  temperature,  and  mixing  ratio)  are  vertically  interpolated  to  the  z(- 
levels  using  a  linear  interpolation  method. 

2.  At  each  z;  level,  the  successive  correction  method  is  performed  using  the  MM5  3D 
data  generated  in  equation  6  as  first-guess  background  data. 

3.  Linear  vertical  interpolation  from  z,  to  zst  is  performed  for  each  WRE-N  parameter 
and  grid  point. 
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3.2  Surface  Data  Analysis 

Upon  completion  of  the  upper-air  analyses  described  in  section  3.1,  a  separate  set  of 
surface  analyses  are  perfonned  for  each  WRE-N  parameter.  Using  the  AFWA  operational 
15-km  resolution  MM5  as  background  data,  objective  analyses  of  surface  temperature, 
dew-point  temperature,  and  horizontal  u  and  v  wind  vector  components  are  made.  The 
following  variant  of  the  successive  correction  method  is  used  in  these  analyses. 

For  a  given  analysis  time,  if  (pa  ( i,j )  is  the  analyzed  value  of  some  meteorological 
parameter  at  grid  point  (/,/)  and  (pt{i,j )  is  the  corresponding  MM5  background  value  for 
that  parameter  (at  that  grid  point),  we  can  write 

m 

<P.  M + 1  -  <PJ  (8) 


where  (p  is  the  observation  at  the  kth  location,  ^  ,  is  the  weight  for  each  observation, 
(0  is  the  value  of  the  background  at  the  observation  point  derived  by  a  bilinear 

/  b,k 

interpolation  method,  and  m  is  the  number  of  observations. 

The  weighting  factor  is  defined  as 

Wk,y  =  ~Y~  (9) 

no¬ 
where  yk  is  the  distance  between  the  kth  observation  point  and  the  grid  point  (/,  /). 

The  method  is  repeated  in  an  iterative  fashion,  so  that  the  background  field  is  updated  by 
the  latest  analysis  after  each  iteration: 

m 

k= 1 

where  (p  (n)  is  the  value  of  the  analysis  at  the  grid  point  after  the  nth  iteration,  and 
(O  (n)  is  its  value  interpolated  for  the  ktb  location  after  the  nth  iteration. 

/  a,k 

In  this  development,  the  iteration  was  repeated  three  times,  for  practical  purposes,  to 
obtain  the  final  value. 

3.3  MARIAH  Similarity  Formulae 

The  equations  and  approach  for  detennining  the  similarity  scaling  constants  for  wind, 
temperature,  and  specific  humidity — referred  to  as  MARIAH  (Rachele  et  ah,  1995) — are 
used  to  combine  surface  and  upper-air  analysis  data.  For  practical  purposes,  it  is  assumed 
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that  the  lowest  level  above  ground  of  the  upper-air  analysis  data  is  modified  by  the 
surface  analysis.  Details  of  this  method  have  been  reviewed  by  Henmi  and  Dumais 
(2004). 


4.  Nowcast  Method  to  Complete  WRE-N 


After  the  analysis  (WRE)  component  of  the  WRE-N  (described  in  section  3),  the 
following  procedures  are  used  to  produce  a  nowcast  out  to  3  h  for  the  same  domain,  the 
final  component  of  the  WRE-N  system. 

For  the  current  WRE  time  (time  =  to),  an  arbitrary  parameter,  such  as  temperature,  dew¬ 
point  temperature,  or  u  and  v  wind  vector  components,  is  expressed  by  A0.  where  A 
represents  the  meteorological  parameter.  For  the  to  +  3  h  nowcast  time  (time  =  t3),  only 
the  3D  MM5  forecast  data  is  available  (no  observations  exist  in  the  future).  These  MM5 
forecast  data  are  then  analyzed  at  time  t3  for  the  domain  of  interest,  using  the  methods 
discussed  in  section  3.  An  arbitrary  MM5  forecast  parameter  at  time  t3  is  expressed  by  A3. 

The  time  rate  of  change  of  parameter  A  between  time  to  and  t3  is  calculated  as 

dX  _  X?,~ Xo  /jj, 

dt~  U-U 


If  no  previous  nowcast  has  been  made  within  the  3-h  period  preceding  time  to.  a  nowcast 
for  the  parameter  A  at  time  t  between  t0  and  t3  is  generated  by 


rtX 

Xr-X^P^it-O- 


(12) 


Here,  [>  is  an  empirical  parameter  to  be  detennined  in  section  6. 


If  a  previously  generated  WRE-N  dataset  that  has  been  generated  within  A t  prior  to  the 
currently  desired  WRE-N  time  (to)  exists,  the  most  recently  generated  WRE-N  fields  will 
be  utilized  in  producing  the  new  nowcast.  To  do  so,  new  variables  XAt  are  defined  for  A  at 
At.  The  time  rate  of  change  of  A  from  At  to  t3  is  calculated  as 


(dX\  (XrXj 
\  dt  i  + 


(13) 


and  the  nowcast  of  A  at  time  t  due  to  XAt  is  calculated  as 

f  dX^ 


XrX^P' 


\dt  j 


1  (t  +  At ) . 


(14) 


If  both  X ,  and  X  va'ucs  are  available,  the  new  nowcast  value  of  A  at  time  t  is 
calculated  as 


10 


x_WrXt  +  *Xr 
wt  +  w' 


(15) 


where 


wt  and  w  ’  are  weighting  factors  defined  as 


W,= 


<-t  o 


and 


1 


W  = 


t  +  A  t 


(16) 


(17) 


In  practice,  A t  is  limited  to  less  than  3  h.  If  previous  nowcast  fields  are  3  h  or  older  than 
the  current  desired  WRE-N  time  to,  they  are  not  used.  An  empirical  parameter  [i  is  used  to 
adjust  the  time  rate  of  change  of  X  in  equations  2  and  4.  In  section  6,  reasonable  values 
for  fi  will  be  investigated. 


5.  Computer  Programs  and  Input  Files 


The  present  WRE-N  is  composed  of  the  following  FORTRAN  77  program  and  sub¬ 
programs: 

1 .  nowcast_main.fi  This  is  the  main  module  that  reads  MM5,  upper-air  radiosonde 
sounding  data,  and  terrain  data.  The  subroutines  3dobj.f  3dobj _dpt.f 
3dobj_ncompf  and  nowcast Jmetsf  are  called  from  within  this  program.  After 
upper-air  and  surface  objective  analyses  of  meteorological  data  are  executed  (and 
merged  using  MARIAH  similarity  theory)  for  times  to  and  A  (discussed  in  section 
3),  nowcast  calculations  are  then  produced  as  discussed  in  section  4.  An  output  file 
short  ranee  fcst  fde  is  produced  upon  completion  of  this  program,  representing 
the  final  4D  WRE-N  grid. 

2.  3dobj.fi  This  subroutine  is  used  when  no  upper-air  radiosonde  sounding  data  are 
available.  The  MM5  forecast  data  are  three-dimensionally  analyzed  onto  grid  points 
of  the  WRE-N  for  times  to  and  1 3. 

3.  3dobj_dpt.fi and  3dobj_ncomp.fi  When  upper-air  radiosonde  sounding  data  are 
available  within  the  WRE-N  time  and  space  domain,  the  MM5  forecast  data  are 
composed  with  the  observed  upper-air  radiosonde  data  (discussed  in  section  3). 
These  programs  are  used  only  for  time  to. 

4.  nowcast_imets.fi:  The  observed  surface  meteorological  data  are  objectively  analyzed 
with  a  successive  correction  method  (discussed  in  section  3)  using  the  data  analyzed 
by  3dobj. f  or  3dobj _dpt.fi/3dobj _ncomp.fi as  background  data,  and  are  subject  to  the 
MARIAH  similarity  formulas. 
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To  use  as  much  observed  data  as  possible  per  WRE-N  cycle,  all  surface  observation  data 
observed  within  ±30  min  of  to  and  all  upper-air  radiosonde  sounding  data  taken  within  2 
hr  of  to,  are  used.  Only  gross  quality  control  measures  are  used  at  this  time,  and  currently 
no  special  weighting  is  used  to  differentiate  between  ingest  data  of  different  ages.  The 
present  scheme  is  designed  to  cycle  at  hourly  intervals,  although  30-min  cycling  is 
possible.  Also,  the  present  configuration  of  the  WRE-N  produces  output  fields  at  30-min 
intervals  out  to  3  h. 

5.1  Examples  of  the  WRE-N 

In  this  section,  a  WRE-N  example  valid  at  0000  Universal  Time  Coordinated  (UTC),  10 
May  2005  for  the  previously  described  Dallas  domain  is  shown.  Figure  5  shows  the 
15-km  resolution  MM5  horizontal  wind  vector  distribution  at  the  10-m  level  for 
0000  UTC,  10  May  2005,  prior  to  interpolation  onto  the  liner-resolution  WRE-N  grid. 


Figure  5.  MM5  grid  horizontal  surface  wind  vector  distribution 
at  10  magi  for  0000  UTC,  10  May  2005. 

Figure  6  is  the  15-km  resolution  MM5  horizontal  wind  vector  distribution  at  the  10-magl 
level,  after  interpolation  onto  the  finer-resolution  1.5 -km  WRE-N  grid  for  the  Dallas 
domain.  In  this  figure,  for  graphical  clarity,  vectors  are  drawn  only  at  every  other  grid 
point.  Observed  surface  wind  vectors  at  the  same  valid  time  are  shown  in  figure  7. 
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Figure  6.  The  10-magl  horizontal  surface  wind  vector  distribution 
interpolated  from  the  MM5  grid  to  the  WRE-N  grid  for 
0000  UTC,  10  May  2005. 


NOTE:  For  clarity,  only  vectors  at  every  other  grid  point 
are  drawn. 


Figure  7.  Observed  10-magl  surface  wind  vectors  valid  at 
0000  UTC,  10  May  2005. 
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Maximum  and  minimum  wind  speeds  obtained  from  just  the  interpolated  15-km 
resolution  MM5  forecast  data  for  this  case  are  6.9  and  1.5  m/s,  respectively.  Figure  8 
shows  the  time  t0  WRE  10-magl  surface  wind  field,  valid  for  0000  UTC,  10  May  2005, 
obtained  after  the  various  successive  correction  analyses  and  similarity  theory 
calculations  (to  merge  the  upper-air  and  surface  analyses)  were  executed.  Maximum  and 
minimum  wind  speeds  of  the  WRE  wind  fields  are  increased  from  those  obtained  solely 
from  the  MM5to  8.9  and  2.3  m/s,  respectively.  Notice  that  the  WRE  wind  directions  are 
also  modified.  Figure  9  shows  the  2  h  (to  +  2  h)  nowcast  10-magl  surface  level  wind  field 
valid  at  0200  UTC,  10  May  2005,  based  on  the  WRE  wind  field  (to)  shown  in  figure  8 
(discussed  in  section  4). 


Figure  8.  The  10-magl  surface  wind  field  produced  by  the 

WRE-N  analysis  module  (0000  UTC,  10  May  2005), 
by  combining  the  fields  shown  in  figures  6  and  7. 
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Figure  9.  The  10  magi  surface  wind  field  valid  for  0200  UTC, 

10  May  2005,  the  2-h  gridded  nowcast  produced  by 
the  WRE-N  extrapolation  module. 

The  2-magl  level  surface  temperature  field  generated  purely  from  the  15-km  MM5 
forecast  data  is  shown  in  figure  10.  For  this  temperature  field,  the  maximum  and 
minimum  temperatures  were  32.8  and  30.3  °C,  respectively.  Figure  1 1  shows  the  2-magl 
level  surface  temperature  field  after  the  execution  of  all  the  analysis  steps  necessary  to 
produce  the  to  WRE.  The  2-magl  level  surface  WRE  temperature  field  has  both  maximum 
and  minimum  temperatures  decreasing  to  27.0  and  22.5  °C,  respectively.  This  is  more 
than  5  °C  lower  than  temperatures  predicted  by  the  operational  AFWA  MM5  15-km 
resolution  forecasts. 
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Figure  10.  The  2-magl  surface  temperature  distribution 

from  the  MM5  interpolated  to  the  WRE-N  grid, 
valid  at  0000  UTC,  10  May  2005. 


Figure  11.  The  2-magl  surface  temperature  field  produced 
by  the  WRE-N  analysis  module  by  combining 
the  temperature  field  shown  in  figure  10  with 
observations  valid  for  00  UTC,  10  May  2005. 

Figure  12  shows  vertical  distributions  (valid  at  0000  UTC,  10  May  2005)  at  the  KFWD 
(Fort  Worth,  TX)  radiosonde  site  for  temperature  and  dew-point  temperature.  Three 
vertical  profiles  are  plotted,  generated  from  the  original  15 -km  resolution  MM5  forecast, 
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the  raw  upper-air  radiosonde  sounding  observation,  and  the  WRE.  The  vertical  profiles  of 
wind  speed  and  wind  direction  for  the  same  location  and  time  are  shown  in  figure  13. 


Figure  12.  Vertical  profiles  of  temperature  and  dew-point  temperature  valid 
for  the  KFWD  (Fort  Worth,  TX)  radiosonde  site  at  0000  UTC,  10 
May  2005.  Profiles  produced  by  the  AFWA  MM5,  the  KFWD 
radiosonde,  and  the  WRE-N  analysis  module  are  shown. 


Figure  13.  Same  as  figure  12,  except  for  wind  speed  and  wind  direction. 
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5.2  Preliminary  Evaluation  of  the  WRE-N  Method 

All  1 5 -km  resolution  MM5  data,  hourly  surface  observation  data,  and  upper-air 
radiosonde  sounding  data  were  collected  over  the  Dallas  domain  shown  in  figure  3  for  the 
period  from  16  April-12  May  2005.  In  this  1 1 1  km-by-1 1 1  km  area  (75-by-75  grid  points 
with  a  1.5 -km  grid  resolution),  there  are  16  regularly  reporting  surface  stations. 

Similar  datasets  were  collected  for  the  period  from  16  May-25  May  2005  over  the 
Denver  domain  shown  in  figure  14,  covering  an  area  of  248  km  by  248  km  ( 125-by- 125 
grid  points  with  a  2-km  grid  resolution)  centered  at  39.5°  N.  and  104.1°  W.  There  are 
only  12  surface  stations  within  this  domain,  all  over  the  plains  area  to  the  east  of  the 
Rocky  Mountain  frontal  range. 

Since  there  are  only  a  limited  number  of  surface  and  upper-air  radiosonde  stations 
available  in  both  domains,  all  the  observations  are  used  for  producing  the  WRE  analyses 
at  time  to  (0  h).  In  the  remainder  of  this  section,  all  available  surface  observations  are  also 
used  to  evaluate  the  quality  of  the  WRE-N  fields  at  times  to,  to+1  h,  t0+ 2  h,  and  /q+3  h. 
Temperature  and  dew-point  are  evaluated  at  the  2-magl  level,  and  winds  are  evaluated  at 
the  10-magl  level.  Statistical  parameters,  such  as  mean  absolute  difference  (AD)  and 
correlation  coefficient  (CC),  are  calculated  at  each  WRE-N  time  period  for  the 
parameters  of  temperature,  dew-point  temperature,  and  wind  speed.  Root  mean  square 
vector  error  (RMSVE)  and  mean  wind  direction  difference  (MWDDF)  are  also  calculated 
for  the  horizontal  wind  vector.  These  statistical  parameters  are  described  in  Efenmi 
(2002). 


Figure  14.  Terrain  contour  distribution  for  the  125-by-125 
domain  of  2-km  grid  spacing,  centered  at 
39.5°  N.  and  104.1°  W.  (Denver  domain). 

NOTE:  The  symbol  represents  MM5  grid  points,  and  the 
“S”  and  “U”  symbols  represent  the  locations  of  the  surface  and 
upper-air  radiosonde  meteorological  observations,  respectively. 
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WRE-N  calculations  with  different  values  of  /?  (0.0,  0.2,  0.4,  0.6,  0.8,  and  1.0),  along 
with  WRE-N  produced  with  no  observation  data  (MM5  data  interpolation),  are  compared 
with  observations  for  the  surface  fields.  A  /E value  of  0.0  represents  persistence.  On  the 
other  hand,  a  /j-valuc  of  1.0  represents  total  reliance  on  the  pure  AFWA  MM5  forecast  at 
the  end  period  (to  +  3  h). 

Figures  15  and  16  show  the  time  series  of  AD  and  CC  for  surface  temperature  for  both 
the  Dallas  and  Denver  domains,  respectively.  For  the  Dallas  domain,  the  current  WRE-N 
method  produced  better  statistics  than  relying  solely  on  the  existing  AFWA  MM5 
forecasts  for  up  to  3  h  forward  from  the  WRE  time  to.  However,  for  the  complex  terrain 
of  the  Denver  domain,  at  the  WRE-N  time  of  t0+3  h,  the  statistics  of  the  AFWA  MM5 
exceeded  those  obtained  by  the  current  WRE-N  method  (for  all  ^-values  tested).  For  the 
Denver  domain,  the  current  WRE-N  method  produced  better  statistics  than  those  from  the 
AFWA  MM5  up  until  to+2  h.  It  should  also  be  noted  that  the  statistics  for  surface 
temperature  produced  from  the  AFWA  MM5  over  the  Denver  domain  are  better  than 
those  produced  over  the  Dallas  domain. 

Figures  17  and  18  are  similar  to  figures  15  and  16,  except  they  show  surface  dew-point 
temperature.  For  the  Dallas  domain,  a  persistence  forecast  (ft  =  0.0)  produced  the  best 
statistics  and  a  /?  =  1.0  produced  the  worst  statistics,  indicating  that  over  this  domain  dew¬ 
point  temperature  changes  slowly  with  time  (and  may  not  be  forecasted  well  by  the 
AFWA  MM5).  On  the  contrary,  for  the  Denver  domain,  the  best  surface  dew-point 
temperature  statistics  were  obtained  with  a/?  =  1.0.  Over  the  Dallas  domain,  the  present 
WRE-N  method  could  be  statistically  better  than  the  pure  AFWA  MM5  forecast  up  to 
to+  3  h,  but  over  the  Denver  domain,  the  same  WRE-N  method  was  statistically  better 
than  the  AFWA  MM5  only  out  to  to+  2  h. 

Figures  19  and  20  display  similar  statistics  for  wind  speed.  It  is  noted  that  the  surface 
wind  speed  statistics  of  the  AFWA  MM5  forecasts  over  the  Denver  domain  are  better 
than  those  produced  over  the  Dallas  domain.  Over  the  Dallas  domain,  the  present 
WRE-N  method,  with  various  values  of  P,  yielded  better  statistics  than  the  AFWA  MM5 
out  to  /o+  3  h,  with  ft  =  0.2  producing  the  best  statistics.  This  may  indicate  that  wind 
speed  changes  slowly  with  time  over  this  domain  and  is  more  uniform  in  distribution. 
Over  the  Denver  domain,  the  method  using  /?  =  1.0  produced  the  best  statistics,  with 
ft  =  0.0  generating  the  worst  statistics  for  wind  speed.  Again,  this  may  be  a  consequence 
of  the  surface  wind  field  varying  significantly  across  the  Denver  domain  due  to  the 
complex  topography,  and  the  AFWA  MM5  model  may  have  some  skill  at 
predicting/resolving  the  local  diurnal  mesoscale  wind  field. 

The  RMS VE  and  MWDDF  surface  wind  field  statistics  are  plotted  in  figures  2 1  and  22 
for  the  Dallas  and  Denver  domains,  respectively.  For  wind  vectors,  the  AFWA  MM5 
forecast  is  considerably  better  for  the  Dallas  domain  than  for  the  Denver  domain.  For 
both  domains,  the  RMSVE  and  MWDDF  statistics  generated  by  the  present  WRE-N 
method  (with  various  ft  values)  are  smaller  than  those  from  just  the  AFWA  MM5 
forecasts  out  to  tu+  3  h.  It  can  be  seen  that  both  the  surface  wind  RMSVE  and  MWDDF 
statistics  for  the  Dallas  domain  are  smaller  than  those  for  the  Denver  domain. 

From  these  results,  it  can  be  tentatively  concluded  that,  for  practical  purposes,  an 
appropriate  value  for  ft  is  0.5.  With  further  study  and  added  complexity,  this  parameter 


19 


could  be  allowed  to  vary  based  on  a  particular  meteorological  and  geographical  scenario, 
though  to  do  so  in  an  automated  fashion,  without  user  input,  would  be  challenging. 


Temperature  (Cl 


Figure  15.  AD  and  CC  of  surface  temperatures  as  a  function 
of  the  WRE-N  nowcast  period  for  different  ft 
values  and  for  the  Dallas  domain. 

NOTE:  The  /I  values  are  given  at  the  right  ends  of  lines. 

The  line  with  “MM5”  refers  to  AFWA  MM5  forecast  data 
(i.e.,  the  background  fields  for  the  WRE-N). 


Figure  16.  Same  as  figure  15,  except  for  the  Denver  domain. 
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Figure  17.  Same  as  figure  15,  except  for  the  dew-point  temperature. 


Figure  18.  Same  as  figure  17,  except  for  the  Denver  domain. 
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Figure  19.  Same  as  figure  15,  except  for  the  wind  speed. 


Wind  Speed  <m/e) 


Figure  20.  Same  as  figure  19,  except  for  the  Denver  domain. 
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Figure  21.  The  RMSVE  and  MWDDF  as  a  function  of  the 
WRE-N  nowcast  period  for  the  Dallas  domain. 


RMSVE  &  MWDD3F 


Figure  22.  Same  as  figure  21,  except  for  the  Denver  domain. 
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8.  Conclusion 


An  initial  WRE-N  method  for  the  IMETS  has  been  developed,  discussed,  and  evaluated. 
In  the  present  method,  meteorological  parameters  including  temperature,  dew-point 
temperature,  and  u  and  v  horizontal  wind  vector  components  (and  surface  pressure, 
although  it  was  not  discussed)  are  analyzed  and  nowcasted.  At  first,  the  3D  fields  are 
computed  across  the  entire  WRE-N  3D  grid  for  each  meteorological  parameter 
(arbitrarily  represented  by  X0  and  Xj)  at  times  to  and  to  +  3  h,  respectively.  The  nowcast  of 
these  parameters,  at  times  between  t0  and  to  +  3  h,  are  done  under  the  assumption  that 
these  parameters  change  linearly  with  the  time  rate  of  change,  multiplied  by  an  empirical 
parameter  [>.  Additionally,  in  cases  where  a  WRE-N  parameter  field  (Xm)  calculated  from 
a  previous  time  At  (less  than  3  h  prior  to  t0)  exists,  the  WRE-N  is  based  on  A0,  A3  and  AA/. 
Only  the  most  recent  WRE-N  field  prior  to  the  t0  is  used  in  the  calculations. 

This  method  has  been  applied  to  two  different  areas  centered  near  Dallas,  TX,  and 
Denver,  CO,  and  has  been  evaluated  using  surface  observation  data  from  the  spring  of 
2005.  It  was  shown  that  an  appropriate  value  for  the  empirical  parameter  [:>  varies 
between  meteorological  parameters  and  different  domains.  An  appropriate  value  for  [i 
was  tentatively  determined  as  0.5,  although  the  statistical  evaluation  studies  described  in 
this  report  are  based  on  the  comparison  between  WRE-N  values  and  observation  for  only 
limited  seasonal  periods  and  geographical  areas.  This  evaluation  study  was  done  only  for 
cases  where  a  previous  WRE-N  field  was  not  used  to  generate  a  new  WRE-N.  Further 
studies  are  needed  to  evaluate  the  present  method  for  cases  in  which  previous  WRE-N 
fields  are  utilized.  The  present  forecast  method  was  shown  valid  up  to  3  h  over  the  Dallas 
domain,  and  up  to  2  h  over  the  Denver  domain.  Again,  tentatively,  we  conclude  the 
method  can  be  valid  for  producing  a  2-  or  3-h  WRE-N,  adding  some  value  over  relying 
on  the  existing  background  AFWA  15-km  MM5  forecast  data  alone. 

Because  there  is  no  physical  and  dynamic  coupling  of  the  WRE  analysis  to  an  NWP 
model  when  producing  the  current  WRE-N  fields,  certain  meteorological  situations  (such 
as  those  with  fast  moving  meso/micro  boundaries  and  thunderstorm  gust  fronts)  will  pose 
significant  difficulty  unless  the  AFWA  MM5  forecast  can  capture  these  features 
accurately  in  time  and  space.  To  address  this  problem  (as  well  as  to  provide  a  capability 
to  ingest  many  other  types  of  non-conventional  meteorological  observations),  a  more 
sophisticated  WRE-N  is  being  envisioned  for  the  Distributed  Common  Ground  Station- 
Army  (DGCS-A)  for  development  later  this  decade.  The  approach  being  developed  is 
based  on  diabatic  initialization  and  would  be  a  fully  coupled  analysis  and  NWP 
assimilation  system.  At  this  time,  a  LAPS  (Shaw  et  ah,  2004)  and  WRF  (Skamarock  et 
ah,  2005)  coupled  system  is  being  investigated  for  application  as  this  advanced  WRE-N. 
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Acronyms 


3D 

three-dimensional 

AFWA 

Air  Force  Weather  Agency 

AOI 

area  of  interest 

ARL 

U.S.  Army  Research  Laboratory 

CC 

Correlation  Coefficient 

DGCS-A 

Distributed  Common  Ground  Station-Army 

DMA 

Defense  Mapping  Agency 

DTED 

Digital  Terrain  Elevation  Data 

GRIB 

gridded  binary 

IMETS 

Integrated  Meteorological  System 

MM5 

Mesoscale  Model  Version  5 

MWDDF 

Mean  Wind  Direction  Difference 

NWP 

Numerical  Weather  Prediction 

RMSVE 

Root  Mean  Square  Vector  Error 

UTC 

Universal  Time  Coordinated 

UTM 

Universal  Transversal  Mercator 

WMO 

World  Meteorological  Organization 

WRE 

Weather  Running  Estimate 

WRE-N 

Weather  Running  Estimate-Nowcast 
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